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Proton/surface interaction pathways: 

• The end state for proton in surface 
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interactions can be one of the 
following [Starukhina, 2001, 2006] : 

- H thermally migrates & diffuse out of 
crystal to form neutral H exosphere 

p+ emission or H 2 emission 

- H gets trapped in crystal vacancy 

- H creates a local ionization event to 
■ immediately react [Poston et al., 

2012] 

- Proton reflects from surface back into 
solar wind [Saito et al., 2008] 

The observation of substantial neutral 
H emission from the lunar surface 
[McComas et al., 2009] and H 2 [Cook 
et al., 2013] suggests implantation 
and outward diffusion is a substantial 
pathway 



H Implantation & Internal Reactions 
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Starukhina 2006 


From Starukhina [2001, 2006, 2012], an 

implanted H will migrate from the surface 
with diffusion time 

t d = h 2 D 0 1 exp(U/kT) 

- Depth: h ~10' 7 m, 

- Diffusion D = D 0 exp (-U/kT) 

- U is diffusion 'activation' energy defining 
both mobility and trapping of H at a position. 

- With defects, U has a range of values, but U 
> 1 eV in regions leading to trapped or 
'loitering' H's in vacancies 





H Diffusion (or Residency) time 
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Farrell et al., 2014 

Diffusion Time Vs temperature 
H Retentive Surfa) 

Vacancy-type Defects 

U = 1.2 eV 



Temperature (K) 

Diffusion time as a function of temperature 


Diffusion (outward) time, 
x D = h 2 D 0 ' 1 exp(U/kT) can 
also be considered the H 
residency time 

Residency time is STRONGLY 
controlled by temperature 
and activation energy, U 

High U: the H can 'loiter' for 
10's of days to possibly form 
OHs..H retentive surface 

Low U value: the H is 
emitted back out of the 
surface in seconds ..H 
emitting surface 


for a family of activation energies, U. 


Number of Implantations per dU 


New Statistical Approach: 
Distribution of Activation Energy 


F(U) ~ exp(-(U-U c )/U w ) 



U c = distribution center, 
U w = distribution width 


For a given sample, the numerous 
implantations have a distribution 
of U values - not a single U but an 
average <U> and 5U deviation 
about the average. 

Expect many of implantations into 
complete lattice structure (with 
low U), but some implantations 
onto lattice with defects 
represented by higher U 

This statistical view lends itself to a 
Monte Carlo approach 

We attempt to search for the 
distribution of U's and 
temperatures that provides long- 
lasting, loitering Hs ....candidates 
to form OH. 


• Divide the lunar equatorial 
or 15° long bins) 

• Solar wind ions flux represented by N = 240000 test particles delivered to 
surface in a flux which varies with local time (i.e., cos of SZA) 

• For each implantation, assign an individual U value quasi-randomly but 
weighted by a pre-defined Gaussian. 

• We then determine diffusion (residency) time as a function of local surface 
temperature. 

• The number of H's with residency time,x D > 10 5 s are considered loitering Hs 
and candidates to eventually form OH. We count these in each local time 
bin. 

- H Continuity Eq.: dn H /dt= S- L diff - L 0H 

- Cold, terminator: Balance near-terminator solar wind flux, F, onto surface vs 
estimated amount of OH content from IR observations: 

T > n 0H h/F > 10 5 s (cannot be much longer than a % lunation) 

• We thus examine the conditions to obtain substantial amounts of these 
'loitering', long-lasting Hs that are candidates to form OH. 

• Example cases now shown 


ethodology 

frontside region into local time bins (~1 hour LT 


Case A: U r = 0 eV , U 
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:a U e = 0.0 eV, U w = 0.1 eV 


0.4^ 


U for Implanted H 
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= 0.1 eV 

Nearly defect free 
lattice with 
low U values 

No loitering H's! 


Implantation Number 

-Most of implantations have activation energies below 0.4 eV. 

-Retention time/loitering time of H's less than 1 second at all locations - all H's lost, no OH. 





Case B: U r = 0 eV , U = 0.3 eV 


U c = 0.0 eV, U w = 0.3 eV 

U for Implanted H 


ing H's 
ors 



a 1.0 x 10 s 2.0 x 10* 8 10 12 TT 

Implantation Number Local Time 

-Only at the terminators is there a substantial population of loitering Hs 

-The warmer regolith near noon releases the Hs, possibly creating a neutral H exosphere 





CaseC: U = 0 eV , U w = 0.6 eV 
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Diffusion Time for Implanted H 

Jo 

. .. . ' : -i 



“O 

o 

c 

ro 
■*— » 
QJ 
Cd 

i_ 

d> 

_Q 


1000- 


Net Retained H vs LT 



No noontime minimum 


0. 1.0 x TO 5 2.0 x 10 s '6 8 10 12 14 

Implantation Number Local Time 

-Substantial population of loitering Hs even at warm noontime 

-Still noontime loss via diffusion, but the loss is offset by greater retention at numerous defects sites 






Comparison to the Moon 
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Clark et al. 2009: Mild retention of 
OH (well below saturation) at 10- 
1000 ppm (over long term) 

Sunshine et al 2009: Reduction in 
OH at sub-solar point 

Consistent with H implantations 
having U c < 0.3 eV, U w < 0.2 eV 


Table 1 

Values of total fraction retained and the relative retention at noon LT (compared to 3 h LT) for a set of runs under varying U t and U w fora thermal 
profile like that at the lunar equator. Bold values indicate UJU W values where a local noontime minimum is modeled. The shaded values are UJU W 
values that could be similar to the Moon's mild retention (see text for further discussion). 



Uc = 0.0 eV 

U c = 0.3 eV 

U ( = 0.6 eV 

Uc = 0.9 eV 

U w = 0.1 eV 

0.0 1 0.0 

0.0 1 0.0 

0.001 1-1*0 

0.22 \-i)M 

U 1V = 0.2 eV 

0.0| 0.0 

0.001 1 - 0.64 

0.034 1 - 0.51 

0.34 |+0.06 

U w = 0.3 cV 

0.001 1 - 0.56 

0,013 [- 0.29 

0.10 | - 0*004 

0.39 1+0.35 

U w = 0.4 cV 

0 . 016 |- 0.19 

0,04 1+0.02 

0,17 |+0.27 

0.42 1+0.51 

U,v= 0.6 eV 

0.10 1+0.35 

0,14 1+0.41 

0.27 1+0.53 

0.45 1+0.67 

Uw=0.75 eV 

0.19 1+0.54 

0.22 |+0.56 

0.33 |+0.63 

0.46 1+0.72 

U w = 0.9 eV 

0.27 1+0.64 

0.30 1+0.65 

0.37 1+0.70 

0.43 1+0.76 



3 Classes of Implantations 


U (eV) Diffusion Loss from Observed Observed 

time (s) surface via in via 


Low U < 0.3 < 10 1 Quick diffusion Emitted H or IBEX, LAMP 

H 2 

Intermediate U 0.3 to 0.9 10 2 -10 14 Diurnal -scale Surface IR(e.g., M 3 ) 

thermal diffusion 

High U >0.9 10 7 To Impactors Surface IR(e.g., M 3 ) 

over 10 16 - 5 


- Model a continuum in U but 'classes' based on loss process 

-We propose that impactors as a loss process regulate the amount of long-term 

loitering H 

-Micro-meteoroids act to 'reset' the affected regolith by releasing H's and by 
adding defects to make the surface more retentive. 

-Keep surface from becoming saturated by 'impact-reset' (Over long term, see 
only 10-1000 ppm -well below saturation) 

-Look at U/T: Explains both thermal effect and increase H in mature (highland) 
soils 








Illustration of Implantation history of a 

1 mm x 1mm area 



Impactor is r ~ 0.1 mm, m = 10' 8 kg Time-> 

Impacts a 1 mm x 1mm region every 10 5 years. 





Conclusions 


• The defect properties of a crystal has as much (or more) 
control on solar wind H retention in an exposed regolith as 
temperature and solar wind influx 

• U/T is the primary variable in residency, not just T 

- T varies by a factor of 4, U varies by a factor of 100 

- T appears as a controlling variable only at intermediate U values 

• The nature of the crystal may be a variable in H retention and 
OH formation that is under-appreciated. 

• Focus on U/T: Can explain a diurnal cycle to H and also 
increased retention in mature soils (lunar highlands) 

• We use Gaussian's to illustrate of effect but real crystal defect 
distribution (U distribution) likely more complicated (e.g., 
amorphized rims - Noble et al., 2005) 

• Starukhina's work pioneering 

• Extended discussion in recent Icarus paper! 
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Interesting Test Case 
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U c = 0.9 eV, U w = 0.1 eV 


Diffusion Time for Implanted H 
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Comparison to the Moon 
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Clark et al. 2009: Mild retention 
(well below saturation) at 10-1000 
ppm (over long term) 

Sunshine et al 2009: Reduction in 
OH at sub-solar point 

Consistent with H implantations 
having U c < 0.3 eV, U w < 0.2 eV 


Table 1 

Values of total fraction retained and the relative retention at noon LT (compared to 3 h LT) for a set of runs under varying U t and U w fora thermal 
profile like that at the lunar equator. Bold values indicate UJU W values where a local noontime minimum is modeled. The shaded values are UJU W 
values that could be similar to the Moon's mild retention (see text for further discussion). 
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• Mild retention (well below saturation) 

• Reduction in H (and hence OH) at sub- 
solar point 

• Consistent with U c < 0.3 eV, U w < 0.2 
eV 

• Have a fraction of implantations with 
0.3 eV < U < 0.9 eV 








Comparison to Moon hs 
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Morning Moon Afternoon Evening 



• Mild retention 

• Reduction in H (and hence OH) at 
sub-solar point 

• Consistent with U c < 0.3 eV, U w < 
0.2 eV 

• Have a fraction of implantations 
with 0.3 eV < U < 0.9 eV 
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